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Current  =  1  to  30kA 
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Diode  impedance  data  and  temporally  resolved  spectroscopic  measurements  both 
indicate  that  anode  plasma  effects  play  an  important  role  in  these 
experiments . 

We  have  constructed  a  neutral  beam-ablation  plasma  experiment  which 
employs  a  Q-switched  ruby  laser  to  independently  generate  a  dense  ablation 
plasma.  The  neutral  beam  is  generated  by  a  (5-20  keV)  duopigatron  ion  source 
with  charge  exchange  neutralization  cell.  Experimental  results  of  beam 
attenuation  show  excellent  agreement  with  predictions  based  on  available 
cross-section  data.  In  studying  the  attenuation  of  the  neutral  beams  we  have 
considered  primarily  charge  exchanee  processes  although  we  are  estimating  the 
relative  contributions  of  momentum  scattering  and  impact  ionization.  Line 
density  values  from  a  self-similar  expansion  model  show  good  agreement  with 
attenuation  data.  Spectroscopic  measurements  of  the  ablation  plasma  yield 
singly  and  doubly  ionized  lines. 

We  have  developed  a  one-dimensional  computational  physics  code  which 
simulates  the  ablation  of  a  carbon  target  by  a  proton  beam.  This  is 
accomplished  by  coupling  a  sophisticated  col  1 isional-radiative  equilibrium 
(CRE)  ionization  dynamics  model,  a  one-dimensional,  single  temperature  and 
fluid  hydrodynamics  model,  and  an  energy  deposition  model.  The  total 
package  is  then  applied  to  simulating  the  beam-target  interaction  under  a 
variety  of  beam  conditions.  Other  applications  considered  include  heating 
due  to  inner-shell  photon  absorption  and  calculation  of  the  amount  of  K 
radiation  emitted  by  the  direct  interaction  of  a  beam  proton  with  a  target- 
atom. 

The  basic  approach  involves  coupling  a  one-dimensional  hydrodynamic 
code  to  an  ionization  dynamics  and  an  energy  deposition  model.  .Because  the' 
radiation  emitted  by  the  plasma  is  an  important  diagnostic  tool,  we  have 
tried  to  model  the  ionization  dynamics ( very  carefully.  This  is 
accomplished  by  developing  an  ionization  dynamics  model  which  considers  a 
large  number  of  excited  states  and  atomic  processes. 

The  atomic  processes  which  are  considered  in  our  model  are: 

1)  collisional  ionization 

2)  three  body  recombination 
collisional  excitation 
collisional  de-excitation 
spontaneous  emission 
radiative  recombination 
dielectric  recombination 

A  model  to  calculate  the  stopping  power  for  energetic,  '"structured" 
(those  carrying  electrons  into  the  target)  low  Z  ions  has  been  developed. 

Low  Z  ions  are  both  energetic  and  structured  at  typical  specific  energies  of 
a  few  hundred  keV/amu  (e.g.  1  MeV  helium).  The  atomic  structure  is  taken 
into  account  by  describing  the  electron  cloud  of  the  projectile  with  elastic 
form  factors  in  the  first  Born  approximation.  However,  at  these  energies 
the  first  Born  approximation  begins  to  break  down.  Hence,  two  corrections 
are  included:  (i)  the  Bloch  correction  to  account  for  the  scattering  of  wave 
packets  (instead  of  plane  waves)  for  close  collisions,  and  (4i)  the  Barkas 
effect  which  corrects  for  deviations  from  pure  Rutherford  scattering  due  to 
the  displacement  of  an  oscillating  electron.  Structure  is  introduced  more 
approximately  in  these  two  terms.  An  attempt  has  been  made  to  allow 
sufficient  generality  to  describe  ionized  targets  as  well  as  neutral  solids. 
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Introduction 


This  final  report  is  divided  into  several  sections.  The 
first  section  covers  the  (5  month)  Research  Option: 
"Electromagnetic  Radiation  Emission  from  Electron  Beam  Induced 
Ablation  Plasmas."  The  remainder  of  the  report  covers  the 
progress  made  during  a  12  month  period  in  research  concerning 
the  interaction  of  neutral  and  ion  beams  with  solids  and  laser 
induced  ablation  plasmas.  Both  theoretical  and  experimental 
results  are  presented. 


R.O.-I  Research  Option:  Electromagnetic  Radiation  Emission  from 

Electron  Induced  Ablation  Plasmas 

During  the  past  year  the  UM  Intense  Energy  Beam  Interaction  Laboratory 
has  developed  a  state-of-the-art  accelerator  for  the  generation  of  intense 
microsecond  pulselength  electron  beams.  This  generator,  the  Michigan 
Electron  Long  Beam  Accelerator  (MELBA),  operates  within  the  following 
range  of  parameters: 

Pulselength  *  0.1  to  4.0  us,  with  fully  adjustable  pulselength  over  the 
first  1.0  us, 

Voltage  *  -1MV  maximum,  flat  to  within  <  ±10%  over  the  first  1.0  us, 

Current  -  10  kA  or  greater. 

Total  delivered  energy  for  a  2.0  us  pulse  >  10  kJ. 

In  order  to  perform  studies  of  electron  beam  pumped  optical  emission, 
we  have  performed  extensive  research  into  the  problems  of  field-emission 
cathode  designs  for  us  electron  beams  [l].  The  goal  of  this  study  was  the 
generation  of  uniform,  reproduceable  electron  beams  for  long  pulse  lengths 
without  diode  shorting  from  cathode  plasma  closure.  The  most  important 
issue  was  the  gap  plasma  closure  velocity  which  determined  the  maximum 
attainable  pulselengths.  In  the  course  of  our  investigation  we  have 
discovered  that  anode  plasmas  can  also  have  a  significant  effect  on  diode 
performance.  The  following  discussion  is  arranged  in  a  loose  chronological 
order. 
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Much  of  the  diode  data  was  obtained  prior  to  our  spectroscopic 
investigations  and  therefore  represented  a  starting  point  for  this  work. 

The  basic  configuration  for  the  initial  diode  experiments  is  depicted  in 
Figure  I.  Both  carbon  fiber  and  velvet  cathodes  were  studied  at  gap 
spacings  of  6,  8,10  cm.  Typical  experimental  data  for  two  shots,  M486  and 
M491,  with  an  AK  gap  spacing  d^  =  10  cm,  and  a  carbon  fiber  brush 

cathode  of  radius  rc  -  7  cm  is  presented  in  Figure  2.  The  data  has  been 

divided  up  into  two  categories;  type  B,  characterized  by  very  fast  rates  of 
impedance  collapse  and  current  increase,  and  type  A  with  slowly  varying 
character.  One  of  the  important  conclusions  of  this  initial  research  was 
that  a  large  area  carbon  brush  cathode  allowed  significantly  longer 
pulselengths  than  the  velvet  cathodes  because  of  lower  cathode  plasma 
expansion  velocities.  Velvet  cathodes  always  gave  type  B  behavior.  Brush 
cathodes  gave  type  A  shots  (M486)  and  shots  with  both  type  A  and  type  B 
character  (M49 1 )  referred  to  as  type  B. 

The  anode  design  also  underwent  a  series  of  modifications.  These  are 
depicted  in  Figure  3.  Configuration  1  was  a  large  carbon  cylinder  with  a  9 
cm  radius.  The  outer  aluminum  surface  was  coated  with  Aerodag  (carbon 
spray)  to  reduce  potential  sputtering  or  else  left  clean.  Pressure  jumps 
after  the  shot  and  large  currents  indicated  that  the  carbon  spray  formed 
ions  in  the  gap  which  enhanced  the  electron  current.  The  appearance  of  arc 
damage  on  the  aluminum  on  the  sharp  edge  around  the  carbon  led  us  to  a 
redesign  of  the  carbon  block  anode  as  we  suspected  that  outgassing  as 
well  as  field  enhancement  under  the  influence  of  the  intense  electron  flux 


"1 

,1 


K 


EXPERIMENTAL  CONFIGURATION  FOR  DIODE  STUDIES 


INSULATORS  OR  CLEAR  RIN6S 
FOR  SPECTROSCOPY 


m 


BRUSH  CATHODE 


APERTURED  ANODE 


I 


m 

mz 


I 

^  APERTURED 
-  BEAM 
URRENT 
MONITOR 


/N 


RETURN  CURRENT  PATH  - ^ 

TOTAL  CURRENT  MONITOR  - ^ 


m 


Figure  1 


NTAL  DATA 

-  7  cm  ,  AK  6AP  -  10  cm 


500  ns/dtv 


VOLTAGE 
(160  kV/dlv) 


CURRENT 
(5.3  kA/dlv) 


HARD  X-RAYS 


APERTURED 
E-BEAM  CURRENT 


Carbon 


Cathode  Stalk 


A:Aerodag  over  AL 
B:  No  Aerodag 


2 


Cathode  Stalk 


A:  Molybdenum 
B:  Carbon 


was  occuring.  In  configuration  2  the  anodes  were  0.3  cm  thick  and  15  cm  in 
radius,  well  removed  from  impinging  electron  flux.  We  tried  molybdenum 
and  carbon  anodes  in  this  configuration.  In  order  to  withstand  the 
multi-kilojoule  electron  beam  pulses,  reduce  the  X-ray  flux,  and  to  reduce 
the  potential  damage  to  the  insulators  and  cathode  in  the  non-crowbarred 
mode  where  arcing  and  subsequent  anode  melting  and  sputtering  may 
occur,  a  solid  poco-graphite  anode  was  employed.  It  is  interesting  to  note 
that  the  diode  current  at  500  ns  with  the  molybdenum  anode  was  6.7  ±  0.3 
kA  compared  to  the  case  for  the  same  size  carbon  anode  at  8.8  +  0.6  kA. 
From  some  of  the  current  modeling  and  energy  deposition  calculations 
following,  it  is  apparent  that  this  can  be  attributed  to  carbon's  greater 
porosity  and  ability  to  absorb  gases  and  impurities  readily.  Ions,  which 
come  from  plasma  formed  at  the  anode  contribute  to  space  charge 
neutralization  at  the  cathode,  allowing  greater  electron  currents.  The 
operating  experience  with  each  anode  design  is  summarized  in  Table  1. 

Some  observations  on  the  experimental  data  in  Figure  2  can  be  made. 
The  rapid  oscillations  in  both  the  apertured  electron  beam  diagnostic  and 
voltage  during  the  time  of  rapid  current  increase  and  voltage  droop 
suggest  some  of  instability  with  counter-streaming  electrons  and  ions  in 
the  diode  [2].  An  effective  gap  closure  velocity  can  be  defined  by  the  total 
time  there  is  a  signal  on  the  apertured  beam  dignostic  (for  a 
noncrowbarred  pulse),  in  this  case  Ueff  =  dAK/tshort  Ueff  Is  proportional  to 

the  rate  of  impedance  collapse  and  is  always  greater  for  a  type  B  shot.  Ueff 

measures  only  the  contribution  to  gap  closure  from  the  cathode  plasma. 
Table  2  summarizes  the  closure  velocity  information  by  type  of  shot  for 
each  cathode  and  gap  spacing.  The  apparent  closure  velocity  Uapp  is 


shots 

aap(cm) 

anode 

1  at  0.5  us(kA) 

384-396 

10 

1 A 

1 0.6  ±  0.6 

397-402 

10 

IB 

9.1  ±  0.5 

403-406 

8 

IB 

14.1  ±2.1 

407-41 1 

6 

IB 

1 7.5  ±  0.8 

412-417 

10 

2A 

6.7  ±  0.3 

418-438 

10 

2B 

8.8  ±  0.6 

439-446 

8 

2B 

12.1  ±  1.2 

447-452 

10 

2B 

7.4  ±0.9 

453-458 

8 

2B 

1 0.9  ±  0.6 

459-466 

10 

2B 

11.4  ±3.9 

(velvet  cathode) 

467-481 

8 

2B 

11.2  ±2.2 

(velvet  cathode) 


Table  1 


discussed  below  in  the  section  on  current  modeling.  Note  that  U^p 

measures  the  total  contribution  to  gap  closure  from  both  the  cathode  and 
anode  plasmas  (if  any).  The  increased  closure  velocity  for  type  B  shots 
seems  to  indicate  that  closure  of  anode  plasmas  across  the  AK  gap  may 
also  be  contributing  to  the  impedance  collapse,  although  increases  in 
current  due  to  other  mechanisms  would  also  be  perceived  as  an  increase  in 
U^.  Our  further  investigations  of  the  diode  physics  concerned  the  issues 

of  current  modeling,  beam  dynamics,  and  the  central  issue  of  ion 
production  mechanisms  and  absorbed  dose  at  the  anode. 


The  physics  of  the  flow  of  charged  particles  in  these  high  voltage 
diodes  is  well  developed.  For  the  current  densities  in  our  diode  the 
Child-Langmuir  model  is  appropriate.  However,  this  model  is 
1 -dimensional  in  nature  and  is  only  strictly  applicable  to  diodes  with 
infinite  (large)  aspect  ratios  (rc/dAK).  For  real  diodes,  particularly  for 

those  in  our  study  (r^/d^*  1 ),  2-d  effects  are  important  and  the  amount  of 

current  from  the  diode  edges  must  be  ascertained.  The  current  model  used 
was  a  combination  of  the  edge  emission  model  of  Parker  et.  al.  [3]  and  the 
technique  of  Poukey  [4]  for  inclusion  of  ion  flows  and  relativistic 
corrections. 

The  total  current  is  given  by  the  sum  of  a  planar  contribution  and  an 
edge  one. 

'total  “  ^planar^epl  +  ^ipP  +  ^edge^eed  +  ^ied^  ^ ^ 


where 


^planar  ”  ^rc2  •  Jepl  ^  Jcl 


jc)  =  2.334  x  10"6V3/2  (3) 

(d0  -  u(t  - 10))2 


Here  jcj  is  the  Child-Langmuir  expression  [5]  for  current  density  modified 


li 


for  gap  closure.  In  that  expression  j  is  in  amps/cm2  for  V  in  volts,  d0  in 
cm.  t0  is  the  diode  opening  time,  d0  is  the  initial  AK  gap  spacing  and  U  is 
an  empirically  determined  closure  velocity.  As  discussed  above  Ueff  is 

appropriate,  and  in  practice  it  is  also  possible  to  vary  U  until  a  best  fit  to 
the  data  is  obtained.  We  refer  to  this  velocity  as  the  apparent  closure 
velocity  Uapp.  These  two  are  not  in  general  equal,  although  when  they  are 

we  can  have  greater  confidence  in  space  charge  models  for  the  flow,  h2  is 
the  correction  factor  due  to  Poukey  that  modifies  jcl  for  the  presence  of 

ions  (leading  to  enhanced  electron  current  density)  and  relativistic  effects 
(leading  to  reductions  in  electron  current  density  *  10%  for  700  kV).  h  is 
given  by 


h 


3 

4 


x/2 

_ _ sittie  Je _ 

[  sin  &  (i  + 13  siV?^m(^/Ycq5  d - 1]^1 

0 


(4) 


In  this  expression  13  is  a  dimensionless  electron  energy,  and  o(  is  a  model 
parameter  varied  to  include  ion  flows  of  variable  magnitude 


j)p  =o(  jeo)[(2me/m,)( 1  +15)1 1/2  ,  G  =  eV/2mec2  (5) 


In  general  the  integral  for  h  must  be  solved  numerically,  in  the  case  where 


I3»0  (nonrelativistic)  and  o(*  0  (no  ions),  we  may  integrate  analytically  to 
find  h*  1  and  we  recover  the  original  results  of  Child  and  Langmuir. 
Langmuir  [6]  treated  the  nonrelativistic  case  with  o(=l  numerically,  a  case 
he  called  bipolar  flow.  Jory  and  Trivelpiece  [7]  solved  the  relativistic  case 
without  ions  exactly  using  a  series  solution.  Each  of  these  limits  are 
contained  in  expression  4  and  have  been  reproduced  in  testing  of  the  model. 

The  edge  current  is  modeled  as  an  angular  portion  of  the  flow  between 
two  concentric  torri,  with  the  cathode  edge  as  the  portion  of  the  inner 
torus  and  the  anode  and  diode  chamber  as  the  outer  torus.  In  the  limit  of 
large  cathode  radius  compared  to  edge  radius  the  flow  between  concentric 
torri  is  approximately  the  same  as  the  flow  between  infinite  concentric 
cylinders,  a  case  where  analytical  series  solutions  have  been  found  by 
Langmuir  and  Blodgett  [8],  We  recover  the  results  of  Parker  for  edge 
emission 


\dge=(0/2n)4ff2r«(do  +  reo).  Jeed  = 


=  2.334  xlO"6V3/2 
(do  +  reo)2A2 


Where  j1b  is  the  Langmuir-Blodgett  result  for  current  density  for  infinite 
concentric  cylinders.  In  this  expression  A  is  given  by 


A  =1f-(2/5)'y2*(ll/120)T3  -(47/3300)'2f4  +  0.00168^5  (8) 


•  y  v  v  ' 

V-  L.V.  k.  .  ’i 


.  -Is.  Lit  U.V  . 


where  '5'  =  lnt(d0  +  reoVrel,  re  =  reo  +  u(t  -  t0)  (9) 


In  expression  9  note  that  the  model  also  takes  into  account  radial 
expansion  of  the  plasma.  This  must  be  modified  for  example  in  modeling  of 
of  diodes  in  a  strong  external  magnetic  field.  The  correction  of  the  edge 
results  for  relativistic  and  ion  effects  by  multiplication  of  by  h2  in 

equation  6  is  only  approximate  here,  whereas  in  equation  2  it  was  exact. 
We  suspect  that  the  corrections  should  be  on  the  same  order  as  for  the 
planar  case,  and  errors  induced  by  not  accounting  for  ions  and  relativistic 
corrections  are  probably  larger  than  those  introduced  by  this  ad-hoc 
correction.  Equation  5  but  with  jeed  and  j)e(j  complete  the  description  for 


edge  emission. 

The  Importance  of  edge  emission  can  be  seen  by  comparing  this  modeled 
edge  current  with  that  from  the  planar  portion  to  obtain  a  scaling 
argument 


*  p  1  anar^  ^  edge  ^p  1  anar^  ep/^edge^eed 


=  20  [(d0  +  reo)rcA2/d02] 

*  rc/do  for  reo  «  do 

In  the  limit  of  large  aspect  ratio  the  contribution  of  edge  emission 
becomes  negligible,  as  we  expect.  For  the  Melba  diodes  rc/d0  *  1  so  edge 

emission  and  planar  emission  contribute  equally. 


A  space  charge  model  can  be  used  to  define  the  apparent  closure 
velocity.  In  principle  one  could  imagine  that  rather  than  use  the 
experimentally  measured  voltage  and  closure  velocity  to  determine  the 
modeled  current,  it  may  be  possible,  given  the  experimental  voltage  and 
current,  to  solve  for  U.  This  procedure  is  not  possible  given  the  full  model 
above,  thus  an  abbreviated  model  for  the  diode  current  was  developed 
which  takes  edge  emission  into  account  by  defining  an  effective  planar 
diode  emitting  area.  The  diode  current  is  given  by 

I  =  ir(Krc)2h22.334x  1  <r6V3/2  ( 1 0) 

Here  Krc  is  the  ‘effective  radius'.  K  is  a  correction  factor  to  take  into 
account  non-ideal  space  charge  behavior.  K  -  1  for  an  infinite  planar  diode. 
Solving  for  d(t)  =  dQ  -  Uapp(t  -  tQ) 

d0  -  U„p(t  - 10)  =  2.708x  t  <T3Krch[  V3/2/l  ] 1  /2  (ID 

Thus  if  the  diode  can  be  modeled  by  space  charge  flow,  a  best  fit  of 
h[v3/2/|]1/2  should  be  a  linearly  decreasing  function  of  time  whose  slope 

yields  Uapp,  the  apparent  closure  velocity.  At  t0  the  following  must  also 

be  satisfied  in  order  that  the  theory  be  consistent 
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d0  =  2708xi0'3Krc{h(V(t0),OO[V(to)3/2/l(to)],/2}fjt  (12) 


which  determines  K.  K  *  1.4  -  1.9  for  the  10  cm  gap,  7cm  radius  brush 
cathode  shots  considered  above.  In  order  to  implement  this  model,  tQ,  the 

diode  opening  time  must  be  specified.  This  can  take  on  many  values 
depending  on  what  definition  of  opening  time  you  subscribe  to.  A 
sensitivity  analysis  shows  that  t0  can  vary  by  over  0.2  ps  and  U  changes 
by  at  most  1 0%. 

The  results  of  applying  the  model  to  shots  M486  and  M491  is  shown  in 
Figures  4  and  5.  Ueff  =  2.6  cm/ps  for  M486  (type  A).  With  this  value  of  U 

and  the  experimental  voltage,  both  Chi Id-Langmuir  models  (solid  line) 
match  the  experimental  current  points  quite  well  at  both  early  and  late 
times.  Incidentally,  =  2.7  cm/ps  for  this  shot,  signifying  that  we 

might  expect  good  agreement.  The  larger  currents  observed  between  0.45 
ps  and  1.2  ps  suggested  the  presence  of  ions  leading  to  bipolar  current 
flows  and  enhancements  of  electron  current.  The  bipolar  current  model  is 
also  shown  and  has  ions  "turned  on"  (c(  nonzero)  at  0.45  ps  and  1  inearl ly 
increasing  until  full  bipolar  flow  (d-  1)  is  reached  at  0.70  ps.  In  fact 
current  increases  above  Child-Langmuir  were  noted  between  0.4  and  0.6  ps 
on  almost  every  shot  taken.  We  felt  this  was  caused  by  exceeding  some 
sort  of  critical  anode  dose  for  plasma  formation.  Since  our  operating 
conditions  were  often  very  close,  even  for  different  gap  spacings,  this 
behavior  might  be  expected  The  important  and  unknown  issue  was  the 
beam  dynamics  which  would  determine  the  beam  current  density  at  the 


•  * 


current  (kA 


anode.  If  the  beam  was  pinching,  high  absorbed  doses  could  be  reached  and 
anode  plasmas  created  easily  giving  significant  plasma  optical  emission. 

In  addition  the  contribution  of  the  impact  ionization  of  the  background 
must  be  evaluated.  This  effect  can  be  significant  for  long  pulselength 
machines  [9], 

The  same  modeling  approach  is  followed  forM491  (type  B)  but  failed  to 
produce  agreement  at  late  times.  Note  that  the  fluctuating  behavior  of  the 
current  model  is  due  to  its  strong  dependence  on  the  voltage  which  is 
drooping  and  very  spikey.  At  early  times,  using  the  effective  closure 
velocity  Ueff  =  3. 1  cm/us,  there  is  both  a  Chlld-Langmuir  plateau  and 

increasing  currents  consistent  with  the  bipolar  model  between  0.38  and 
0.87  ps.  The  qualitative  agreement  at  late  times  with  the  apparent  closure 
velocities  -  5  cm/ps  may  suggest  the  influence  of  anode  plasma 

closure.  However  the  lack  of  complete  agreement  suggests  that  other 
phenomena  need  to  be  considered  Vacuum  arc  behavior  is  ruled  out  on  the 
basis  of  the  constant  impedance  behavior  expected  of  arcs.  The  currents 
are  also  not  consistent  with  other  models  for  high  current  pinching  diodes 
[10],  It  is  evident  from  observing  the  same  type  of  current  behavior  for 
molybdenum  anodes,  that  the  cathode  plasma  also  plays  a  role  in  the 
formation  of  these  very  large  currents.  There  are  extensive  references  in 
the  Russian  literature  [1 1,12]  to  cathode  plasma  instabilities  in  long  pulse 
beams  leading  to  enhancements  over  Child-Langmuir  current  densities. 

The  mechanisim  for  these  instabilities  is  the  depletion  of  the  cathode 
plasma  by  drawing  electrons  off  faster  than  they  are  being  resupplied  by 
plasma  formation.  This  is  critical  for  long  pulselengths  where  the  plasma 
may  not  be  continually  regenerated  when  the  electric  field  and  electric 
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field  rates  are  reduced  after  the  fast  initial  rise.  This  occurs  when  the  gap 
impedance  has  dropped  below  some  critical  value,  depending  of  course  on 
the  cathode  plasma  properties.  Dropping  impedance  can  occur,  of  course, 
because  of  gap  closure,  but  also  in  this  case  after  ions  form  and  enter  the 
gap.  The  electrons  are  drawn  off  at  the  rate  determined  by  the  diode 
impedance,  however,  and  since  the  plasma  is  not  being  resupplied  ,  it  is 
possible  that  gap  closure  may  be  checked  or  slowed.  A  charge  imbalence 
builds  up  in  the  plasma  and  a  positive  ion  sheath  forms  at  the  cathode 
surface  to  compensate,  leading  to  large  electric  fields.  The  electrons  are 
now  accelerated  and  reach  the  cathode  plasma  emission  surface  with  a 
finite  velocity.  The  configuration  is  now  similar  to  a  triode  with  the 
cathode  flare  acting  in  the  role  of  the  grid  in  a  triode.  The  large  electric 
fields  created  at  the  cathode  by  the  ion  sheath  results  in  a  resupply  of 
material  to  the  bulk  of  the  cathode  plasma.  The  net  charge  imbalence 
disappears  and  the  diode  returns  to  its  normal  operating  characteristics. 
This  instability  may  cycle  once  or  more.  These  large  enhancments  over 
even  bipolar  currents  are  interesting.  With  such  rapidly  falling  voltages 
and  increasing  currents,  the  anode  surface  dosage  increases  substantially. 
Further  investigation  is  necessary. 

Beam  Dynamics  for  Anode  Plasma  Studies 


Computer  simulations  and  modeling  of  the  beam  dynamics  were 
completed  using  the  Stanford  Electron  Optics  Code  (SEOC).  From 
experimental  measurements  of  the  anode  current  density  distribution 
during  the  first  0.5  us  we  had  reason  to  suspect  that  the  center  of  the 
beam  was  pinching  slightly  while  the  outer  portions  were  diverging  (see 


Figure  6).  The  code  simulations  agree  with  this  result.  In  order  to  model 
measured  currents  during  the  first  0.5  ps  (before  the  suspected  presence 
of  ions)  it  was  necessary  to  account  for  edge  emission,  just  as  in  the 
analytical  model.  In  this  case  planar  emission  contributed  more  than  edge 
emission  while  for  the  analytical  model  edge  emission  was  slightly 
greater  than  planar  emission.  The  most  important  information  obtained 
from  the  simulation  was  the  current  density  of  the  beam  on  the  anode, 
which  would  allow  us  to  calculate  the  absorbed  dose  and  to  address  the 
machanisim  for  anode  plasma  formation.  For  the  10  cm  gap,  7  cm  radius 
cathode  shown  in  Figures  7  and  8  this  is  20  A/cm2.  The  simulations  with 
SEOC  converge  with  a  great  deal  of  effort  and  are  at  their  limit  of 
usefulness.  The  code  was  not  designed  for  intense  electron  beams  with 
large  self  magnetic  fields.  As  we  will  progress  to  significantly  higher 
current  densities,  it  is  clear  that  another  simulation  code  will  be 
necessary  We  have  obtained  TIGER,  a  Monte  Carlo  energy  deposition  code. 
This  code  coupled  with  an  adequate  beam  dynamics  code  will  allow  us  to 
calculate  absorbed  doses  in  solids,  gases,  and  plasmas  in  many  different 
experimental  situations. 

Calculations  of  spatially  and  time  varying  absorbed  energy  dose  using 
the  model  of  Spencer  [13]  failed  to  produce  the  doses  required  for  plasma 
formation.  Many  authors  have  reached  the  conclusion  that  200  -  700  J/gr 
on  carbon  or  aluminum  will  cause  the  release  of  desorbed  gases  [3,14,15] 
These  gases  are  subsequently  impact  ionized  or  avalanche  ionized  at 
higher  surface  densities.  These  studies  were  all  done  for  short  (~  100  ns) 
high  current  density  (  *  I  kA/cm2)  pulses  The  absorbed  doses  for 
vaporization  of  carbon  to  form  neutrals  or  ions  is  even  higher  (  50  -  180 
kJ/gr)  [16-17], 
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Using  the  current  density  incident  on  the  anode  given  by  SEOC,  the 
average  energy  during  the  first  0.5  us  (700  kV),  and  0.31 1  gr/cm2,  the 
electron  range  for  that  energy,  the  approximate  anode  dose  is 
45t(ps)  J/gr.  At  t=0.45  ps  the  dose  is  20  J/gr  which  is  an  order  of 
magnitude  below  the  minimum  required  dose.  A  simple  calculation  shows 
that  impact  ionization  cannot  contribute  sufficient  numbers  of  ions  to 
allow  bipolar  flow  at  the  diode  pressures  we  operate  at.  It  thus  seems 
that  the  presence  of  ions  is  not  admitted. 

The  experiments  of  Swain  et.  al.  (14]  were  done  with  in-situ  cleaning  of 
aluminum  anodes.  They  note  that  the  doses  for  anode  plasma  formation  for 
anodes  contaminated  by  pump  oil  or  vacuum  grease  are  reduced  by  a  factor 
of  10  or  more.  Poco-graphite  has  about  a  20%  porosity  and  we  suspect  that 
it  is  even  easier  to  bring  gases  or  impurities  out.  We  suspect  that  pump  oil 
and  water  vapor  are  the  largest  contributors  to  the  impurities  in  our 
system.  The  presence  of  pump  oil  in  our  system  is  consistent  with  our 
experimental  procedure  of  roughing  on  the  chamber  as  continuously  as 
possible.  In  a  turbo  pumped  system,  back  diffusion  of  oil  from  the  roughing 
pumps  will  occur  if  the  blades  are  not  spinning  while  in  a  diffusion 
pumped  system  the  diffusion  oil  also  contributes.  One  final  point  is  that 
the  critical  doses  necessary  may  be  reduced  for  longer  pulselengths. 
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The  preliminary  experimental  configuration  for  measuring  spectra  is 
shown  in  Figure  1.  We  used  polished  lucite  rings,  as  shown,  which  allowed 
light  to  pass  through  the  optics  and  back  to  the  screen  room  to  be  focused 
on  the  OMA  slit  (Optical  Multichannel  Analyzer).  This  was  a  configuration 
which  would  allow  us  to  take  spectra  with  as  little  delay  and  change  to 
the  experimental  configuration  as  possible.  There  were  several  drawbacks 
to  these  rings,  however,  the  first  being  that  the  curved  surface  of  the 
lucite  acted  like  a  dispersing  lens  and  the  thickness  as  an  attenuator. 
Secondly  there  is  a  cutoff  of  blue  light  at  about  380  nm  below  which 
nothing  is  transmitted.  Arcing  to  the  lucite  ocurred  during  a  few  shots  and 
dendrites  under  the  surface,  discovered  by  inspection  afterwards  often 
left  one  with  an  unavoidably  ambiguous  interpretation  of  the  spectra 
obtained.  In  addition,  the  copious  edge  emission  led  to  a  large  electron  flux 
near  the  lucite,  although  no  surface  roughness  was  created.  In  addition  to 
discoloration  of  the  lucite  and  possible  formation  of  a  lucite  surface 
plasma  (although  in  fact  no  oxygen  lines  were  observed),  this  may  have  led 
to  an  intense  continuum-like  emission  or  Cerenkov  radiation  when  the  OMA 
was  gated  during  the  pulse.  It  may  in  fact  be  continuum  emission  or 
possibly  bremsstrahlung.  This  feature  did  not  occur  all  the  time  and  we 
were  able  to  obtain  some  spectra  during  the  pulse.  Even  under  these 
restrictive  conditions  it  was  possible  to  obtain  information  which  gave  a 
better  understanding  of  the  processes  that  were  ocurring.  The  spectra 
obtained  required  careful  interpretation. 

Spectroscopy  on  the  diode  chamber  was  made  possible  with  highly 
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polished,  large  diameter  lucite  rings.  Although  one  set  of  lucite  rings  with 
flat  transmitting  surfaces  was  used,  the  spectra  using  curved  surfaces 
were  not  significantly  different.  Due  to  the  dangers  of  misinterpretation 
of  spectra  in  the  analysis,  we  weighted  the  results,  favoring  spectra 
obtained  early  in  the  run  before  significant  lucite  damage  could  have 
occurred. 

The  first  spectra  were  obtained  using  OMA  gates  the  full  width  of 
noncrowbarred  pulses.  With  these  gates,  light  from  a  whole  range  of 
plasma  producing  events  in  the  diode  would  have  been  observed  including 
cathode  and  anode  plasma  light,  plasma  formed  by  the  late  time  pulses  due 
to  the  ringing  stage  [1]  or  the  plasma  formed  from  any  late  time  arcs 
which  may  have  occurred.  Initially  these  spectra  were  hard  to  interpret 
and  no  lines  were  readily  identifiable.  This  was  due  to  the  large  Stark 
shifts  caused  by  the  diode  electric  fields.  In  fact  the  shift  in  the  Hfi  line  in 

Figure  9  and  10  is  about  1 7  X,  which  corresponds  to  a  diode  electric  field 
of  75  kV/cm.  This  is  in  good  agreement  with  the  theoretically  obtained 
electric  fields  of  100  kV/cm  required  to  see  this  shift.  Shifts  were  also 
observed  in  the  C2  spectrum  and  in  the  very  intense  C+-4265  a  line.  No 

neutral  carbon  lines  were  observed,  apparently  because  oi  their  lower 
intensities  and  since  important  lines  may  have  been  hidden  by  the  intense 
C2  structures.  Since  the  spectra  are  shifted,  they  must  have  been  obtained 

during  the  applied  field  t.e.  during  the  electron  beam  pulse 

For  the  pulses  in  Figures  9  and  10,  the  time  of  the  applied  field  is 
shorter  than  is  needed  for  the  cathode  plasma  to  cross  the  gap,  thus  we 
are  observing  anode  plasma  components.  There  are  some  other  interesting 
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features  of  those  two  spectra.  First,  the  cutoff  in  transmission  below 
about  380  nm  is  identified.  A  broad  spectrum  white  light  source  was  used 
to  verify  this  cuttoff.  Next  note  the  increase  in  continuum-like  structure 
from  the  first  (M551)  to  the  second  (M552)  as  it  obscures  details  and 
hinders  resolution.  By  the  sixth  shot  of  that  day  line  emission  was 
completely  obscured.  Finally,  comparing  the  two  spectra  reveals  that  the 
slightly  longer  pulse  length  (M552),  extending  just  into  the  region  of 
current  ramp-up  resulted  in  3x  more  emission  giving  us  evidence  for  the 
evolution  of  desorbed  gases  from  the  anode  at  Vs  0.5  -  0.6  ps.  This  was 
also  confirmed  on  a  different  day  using  short  OliA  gates  immediately 
following  short  crowbarred  pulses.  Shorter  gates  can  be  use  to  exclude 
all  light  except  that  produced  by  the  anode  ablation  plasma.  No  line  spectra 
were  obtained  for  pulses  i  0.4  ps.  The  spectra  obtained  for  pulses  *  0.45  - 
0.7  ps  were  completely  similar  to  the  structure  in  the  presented  spectra 
but  at  relative  amplitudes  i  100.  One  spectra  was  obtained  when  the  OMA 
gate  bracketed  the  late  time,  low  energy,  high  current  pulse  due  to  the  LRC 
ringing  of  the  reverse  charged  stage.  Again  the  same  structure  was 
observed,  indicating  the  anode  as  source,  but  at  a  relative  amplitude  of 
2500,  indicating  that  the  low  energy  optimizes  the  front  surface  dose  by 
minimizing  the  electron  range. 

The  spectra  obtained  from  a  noncrowbarred  pulse  is  given  in  Figure  1 1. 
The  gate  was  moved  forward  to  about  0.8  ps  to  avoid  the  continuum-like 
structure  that  would  obscure  the  lines,  and  indeed  note  that  the 
background  is  flat  The  CH  peak  is  absent,  the  C+  peak  is  more  clearly 
visible,  and  the  HB  peak  is  reduced  in  amplitude  In  fact  on  consecutive 
shots  that  day  the  hydrogen  peak  was  constantly  reduced  in  amplitude  The 


Stark  shifts  obtained  here  are  smaller,  implying  smaller  average  fields. 

A  large  portion  of  the  spectra  in  all  three  cases  remains  unidentified 
due  to  the  difficulty  when  Stark  shifts  are  present  but  also  due  to  the 
resolution  limits  of  the  instrument.  Theoretically  densities  can  be 
estimated  from  broadening  and  tempertatures  from  line  ratios.  These 
possibilities  will  be  investigated  for  this  plasma  but  resolution  limits 
again  make  the  task  difficult.  Continued  investigations  on  the  plasmas 
formed  in  the  diode  will  be  made  using  a  new  aluminum  ring  with  a 
shielded  observation  port.  This  new  structure  will  eliminate  the  ambiguity 
of  the  previous  scheme  and  allow  us  to  optimize  the  light  output  by  using 
thin  quartz  windows.  If  Cerenkov  or  continuum-like  emission  is  still  a 
problem,  a  smaller  area  cathode  can  be  used  to  minimize  the  edge 
emission.  In  addition  we  will  be  able  to  look  right  on  the  surface  of  a 
carbon  target.  Other  experimental  results  at  this  lab  stress  the 
importance  of  looking  right  along  the  surface  at  the  most  dense  part  of  the 
ablation  plasma. 

Conclusions  on  Preliminary  Results  During  Research  Ootion 


We  understand  the  diode  physics  fairly  well  and  further  investigation 
and  optimization  of  plasma  production  within  the  diode  may  help  fill  the 
remaining  gaps.  We  have  formulated  a  tentative  hypothesis  about  the 
plasma  formation  process  which  will  require  some  further  testing 

I )  Contaminated  anode  surfaces  contribute  to  the  space  charge 
neutralization  of  the  beam  on  timescales  ~  0.5  ps. 


2)  Space  charge  neutralization  ieads  to  pinching  to  a  much  greater 
degree  than  in  the  simulations,  increasing  the  local  current 
density  at  the  anode.  This  may  also  lead  to  cathode  plasma 
instabilities  and  enhancements  over  even  bipolar  currents. 

3)  More  plasma  is  formed  from  the  ususal  anode  desorption 
mechanisims  at  longer  times  *  l  -2  ps  For  noncrowbarred 
pulses  the  energy  deposition  approaches  15  kJ,  front  surface 
doses  are  maximized  and  production  of  carbon  by  vaporization 
mechanisims  may  come  into  play. 

Definite  conclusions  about  the  density  and  temperature  of  these 
plasmas  are  not  possible  from  these  preliminary  experimental  results.  The 
measurements  were  all  taken  at  the  maximum  OMA  gain  because  of  the 
small  intensities  detected.  This  was  due  for  the  most  part  due  to  the 
diffusing  and  attenuating  nature  of  the  lucite  ring.  A  preliminary  estimate 
of  the  lower  bound  on  plasma  density  from  space  charge  neutralization 
considerations  of  the  electron  beam  is  3  x  109  cm"3  These  important 
questions  will  be  answered  with  the  greatly  improved  viewing  method  we 
are  pursuing. 

The  most  important  conclusion  we  have  reached  is  that  methods  for 
increasing  the  beam  current  density  at  the  anode  or  carbon  target  must 
certainly  be  investigated.  This  is  required  both  to  form  the  plasma  and  to 
put  it  in  a  highly  ionized  nonequilibrium  state.  In  order  to  maximize  the 
current  density  one  or  more  of  the  following  steps  can  be  taken. 


1  )We  are  in  the  process  of  acquiring  external  field  coils  for  magnetic 
compression  as  well  as  for  extracting  rotating  beams  which  have 
increased  deposition  due  to  their  shorter  penetration  depths. 

2)  From  the  Child-Langmuir  scaling  it  is  evident  that  decreasing  the 
gap  spacing  is  also  an  alternative.  If  it  is  decreased  sufficiently 
it  may  be  possible  to  create  a  seif  pinched  electron  diode, 
early  in  the  pulse.  Pinching  increases  the  current  density  from  local 
area  decrease  but  also  from  shorter  penetration  depths  due  to 
nonperpendicular  incident  angles  We  must  assess  whether  this 
drives  the  generator  harder  than  it  was  designed  for.  At  present  we 
can  probably  reduce  the  gap  by  a  factor  of  2,  increasing  the  current 
density  by  a  factor  of  4. 
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I.  Theoretical  Progress 

la.  Stopping  Power  of  Low-Z  "Structured”  Pro jecti le 
Ions 

A  model  to  calculate  the  stopping  power  for  energetic, 
"structured"  (those  carrying  electrons  into  the  target)  low 
Z  ions  has  been  developed.  Low  Z  ions  are  both  energetic 
and  structured  at  typical  specific  energies  of  a  few  hundred 
keV/amu  (e.g.  1  MeV  helium).  The  atomic  structure  is  taken 
into  account  by  describing  the  electron  cloud  of  the 
projectile  with  elastic  form  factors  in  the  first  Born 
approximation.  However,  at  these  energies  the  first  Born 
approximation  begins  to  break  down.  Hence,  two  corrections 
are  included:  (i)  the  Bloch  correction  to  account  for  the 
scattering  of  wave  packets  (instead  of  plane  waves)  for 
close  collisions,  and  (ii)  the  Barkas  effect  which  corrects 
for  deviations  from  pure  Rutherford  scattering  due  to  the 
displacement  of  an  oscillating  electron.  Structure  is 
introduced  more  approximately  in  these  two  terms.  An 
attempt  has  been  made  to  allow  sufficient  generality  to 
describe  ionized  targets  as  well  as  neutral  solids. 

The  first  Born  calculation  must  account  for  "shell 
corrections"  at  these  energies.  Shell  corrections  are 
energy  dependent  alterations  to  the  stopping  power  due  to 
the  less  effective  contributions  of  the  inner  shell  target 
electrons  as  the  projectile  slows  down.  To  obtain  realistic 
values  of  shell  corrections,  we  follow  Bonderup's1  version 
of  the  local  plasma  approximation  (LPA) .  The  LPA  is  an 


average  of  the  stopping  power  of  a  uniform  density  plasma 


37 


over  the  target  electron  cloud  distribution. 

Bonderup's  treatment  of  the  LPA  utilizes  large  and 

small  velocity  expansions  (projectile  velocity  relative  to 

characteristic  target  electron  speeds) .  The  dominant 

contribution  from  the  large  velocity  expansion  is  then 

compared  to  a  detailed  atomic  calculation  to  introduce 

corrections  to  the  LPA  for  not  describing  the  binding  of  the 

target  electrons.  For  point  particles  (i.e.  no  structure), 

2 

the  atomic  calculation  was  the  Bethe  formula  with  shell 

3 

corrections  calculated  by  Fano  and  Turner. 

The  stopping  power  calculation  for  structured 

projectiles  of  a  uniform  density  plasma  was  indicated  by 

4 

Ferrell  and  Ritchie.  Their  results  are  used  for  the  slow 
velocity  term.  We  performed  the  high  velocity  expansion 

5 

using  Lindhard's  dielectric  function.  The  first  few  terms 
(without  shell  corrections)  of  the  atomic  calculation  for 
structured  particles  were  obtained  by  Kim  and  Cheng. ^  We 
extended  the  work  of  Fano  and  Turner  to  include  the  shell 
corrections.  Since  we  are  generally  interested  in  only  a 
few  bound  projectile  electrons,  we  shall  ignore  projectile 
excitation. 

Comparing  the  two  high  velocity  expressions,  we 
recognize  the  need  to  introduce  the  Lindhard-Scharf f  binding 
parameter  (LSBP)7  in  the  LPA  definition  of  the  Bethe  mean 
excitation  energy.  This  arises  from  Lindhard's  argument  of 
replacing  the  plasma  frequency  with  an  effective  frequency 
given  by  the  LSBP  times  the  plasma  frequency.  Based  on  a 
Thomas-Fermi - 1  ike  model  of  the  target  atom,  Lindhard 


determined  that  the  LSBP  was  a  constant  of  order  2. 

However,  recent  investigations  of  ionized  targets  by  J. 

g 

Peek  have  indicated  that  the  LSBP  must  be  dependent  on  the 
degree  of  target  ionization  such  that 
LSBN(Nt  =  1)  =  LSBP (N^  =  Zfc) 
where  and  Z  are  the  number  of  target  electrons  and 
protons,  respectively.  This  has  raised  questions  about 
using  the  LPA  for  ionized  targets. 

However,  it  is  not  surprising  that  Lindhard's  arguments 
fail  for  one-electron  target  ions  since  it  is  based  on  a 
Thomas-Fermi-like  (many-electron)  model.  Hence,  we  repeated 
his  arguments  with  a  one-electron  model  of  the  atom.  For 
large  Zfc  we  obtain 

LSBN(Nfc  =  1)  2  1.36 

which  essentially  agrees  with  Lindhard's  neutral  LSBP  value 
and  Peek's  scaling.  For  low  Zt  we  should  not  use  Lindhard's 
many-electron  value.  However,  we  may  use  Peek's  scaling  to 
predict  neutral  values.  These  are  plotted  in  Figure  1. 

The  long  dashed  curve  is  Lindhard's  value;  the 

solid  curve  is  our  calculation.  The  short  dashes  with 
asterisks  are  average  values  based  on  experimental  data 

9 

taken  from  Ziegler.  Our  curve  reflects  the  proper 
general  trend  with  the  target  atomic  number.  This  gives  us 
confidence  in  our  values  for  LSBP(Nt  =  1)  and  helps  put  the 
LPA  for  ionized  targets  on  more  firm  ground.  Data  for  gases 
are  shown  since  the  calculations  are  based  on  an  isolated- 
atom  picture  of  the  target  (appropriate  for  gases  or  highly 


ionized  solids).  The  differences  in  magnitudes  are  not 
important  since  the  LSBP  for  neutral  targets  is  determined 
to  give  the  accepted  value  of  the  Bethe  mean  excitation 
parameter . 

For  partially  ionized  targets,  we  repeat  a  qualitative 
scaling  argument  of  Lindhard's7  for  a  target  with  N 
electrons.  This  suggests  the  scaling 
LSBP  (Nt)  =  LSBP  (Nt  =  Zfc) 

The  results  for  =  13  and  81  are  shown  in  Figure  2  along 
with  detailed  atomic  calculations  of  J.  Peek^  (boxes  and 
asterisks).  Curves  1  and  2  are  based  on  an  empirical  fit  to 

O 

Z  =  13  calculation.  Curves  3  and  4  are  our  calculations 
based  on  theory.  Lack  of  agreement  is  not  surprising  given 
the  crude  arguments  used  to  obtain  our  scaling.  It  seems 
that  our  values  should  be  modified  (increased)  to  obtain 
more  realistic  values. 

Term  by  term,  the  high  velocity  LPA  expression  is 
compared  to  the  corresponding  atomic  calculation.  In  each 
case,  replacing  the  plasma  frequency  with  the  LSBP  times  the 
plasma  frequency  is  sufficient  to  give  the  proper  scaling 
with  target  ionization.  The  actual  numerical  magnitude  may 
have  to  be  adjusted  in  each  case  to  correct  for  deficiencies 
in  the  LPA  (only  one  term  in  the  shell  correction  had  to  be 
significantly  altered).  The  close  comparison  with  the 
detailed  atomic  case  gives  us  confidence  in  our  use  of  the 
LPA. 

Structure  of  the  projectile  is  introduced  in  the  Bloch 
correction  using  the  derivation  attributed  to  Lindhard11 
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log  LSBP(N)/LSBP(N=Z) 


based  on  partial  wave  expansions.  A  point-like  effective 

charge  and  effective  Bethe  mean  excitation  parameter  are 

chosen  similar  to  that  suggested  by  Kim  and  Cheng. ®  These 

reproduce  the  high  velocity  atomic  calculation  without  shell 

corrections  (shell  corrections  are  not  treated  in  the  Bloch 

or  Barkas  terms) .  The  resulting  Bloch  correction  agrees  up 

to  a  factor  of  order  unity  inside  a  logarithm  to  a 

1 2 

calculation  based  on  Bloch's  original  derivation  in 

appropriate  limits.  In  this  latter  model,  an  effective 

charge  equal  to  the  projectile  nuclear  charge  Z  is  used  for 

P 

close  collisions  while  the  net  charge  (Z  -  N  )  is  used  for 

P  P 

distant  collisions.  Target  ionization  effects  follow  from 
the  use  of  an  effective  Bethe  mean  excitation  parameter 
which  tends  to  reduce  the  magnitude  of  the  correction  as  the 
target  is  ionized. 

Structure  is  introduced  in  the  Barkas  effect  again  by 

point-like  effective  charges:  Z^  for  close  encounters  and 

(Z  -  N  )  for  distant  ones.  However,  there  is  a  debate  over 
P  P 

whether  any  contribution  results  from  close  collisions.1^ 

Due  to  the  difference  between  the  net  and  nuclear  charges 

(and  how  that  may  be  magnified  by  the  cubic  charge 

dependence  of  the  Barkas  effect),  structured  projectile 

stopping  power  provides  a  good  test  for  the  correct  model. 

Our  calculations  have  indicated  that  it  is  necessary  to 

include  the  close  collisions.  Our  model  follows  that  of 
14 

Deutsch  with  the  impact  parameter  modified  in  the  usual 
semi -classica 1  way  at  lower  energies.  Target  ionization  is 
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taken  into  account  naturally  through  the  LPA  and  the  LSBP. 
We  use  the  same  target  wave  functions  as  used  in  our  LPA  to 
the  first  Born  approximation. 


Figure  3  shows  our  initial  results  for  structured 

projectile  stopping  in  neutral  solids.  In  particular,  the 

results  are  for  1  MeV  helium  in  various  solids.  Our 

1 5 

results  are  given  by  asterisks,  Ziegler's  tabulations  by 
boxes.  The  agreement  is  promising  with  a  maximum  difference 
of  about  16%.  Solid  state  wave  functions  of  Moruzzi,  Janak 
and  Williams^  were  used  to  describe  the  targets.  The 
projectile  mean  equilibrium  charge  was  taken  from 
calculations  by  Ohtsuki.'*'^  A  two  charge  state  approximation 
was  used  for  helium  (He+^  and  He+^)  with  hydrogen-like  wave 
functions  for  He+1.  The  stopping  number  G  is  proportional  to 


with  j^atom  =  target  atom  density 

=  number  of  bound  target  electrons 
e  =  electron  charge 

m  =  electron  mass 

v  =  projectile  velocity 

As  the  target  becomes  ionized,  free  electrons 
contribute  to  the  stopping  power.  This  may  be  calculated 
from  our  high-velocity  expansion  for  bound  electrons  by 
replacing  the  Fermi  distribution  for  a  degenerate  gas  with  a 
Fermi-Dirac  distribution  of  a  nondegenerate  gas  with  finite 
temperature.  This  was  possible  by  applying  the  expansion 
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1 8 

technique  of  Sigmund  and  Fu  to  structured 

projectiles.  For  the  special  case  of  bare  ions,  our  result 
agrees  with  a  temperature  dependent  calculation  of 
Deutsch1^ . 
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I  b. COMPUTER  SIMULATION  OF  ION  BEAM  TARGET  DYNAMICS 
AND  RADIATION  EMISSION 

We  have  developed  a  one  dimensional  computational  physics  code 
which  simulates  the  ablation  of  a  carbon  target  by  a  proton  beam.  This 
is  accomplished  by  coupling  a  sophisticated  collisional -radiative 
equilibrium  (CRE)  ionization  dynamics  model,  a  one  dimensional,  single 
temperature  and  fluid  hydrodynamics  model,  and  an  energy  deposition 
model.  The  total  package  is  then  applied  to  simulating  the  beam-target 
interaction  under  a  variety  of  initial  beam  conditions.  Other 
applications  considered  include  heating  due  to  inner-shell  photon 
absorption  and  calculation  of  the  amount  of  radiation  emitted  by  the 
direct  interaction  of  a  beam  proton  with  a  target  atom. 

The  basic  approach  involves  coupling  a  one  dimensional  hydrodynamic 
codel  to  an  ionization  dynamics  and  an  energy  deposition  model.  Because 
the  radiation  emitted  by  the  plasma  is  an  important  diagnostic  tool,  we 
have  tried  to  model  the  ionization  dynamics  very  carefully.  This  is 
accomplished  by  developing  an  ionization  dynamics  model  which  considers 
a  large  number  of  excited  states  and  atomic  processes. 

The  atomic  processes  which  are  considered  in  our  model  are: 

1)  collisional  ionization 

2)  three  body  recombination 

3)  collisional  excitation 

4)  collisional  de-excitation 

5)  spontaneous  emission 
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6)  radiative  recombination 

7)  dielectric  recombination 

These  are  incorporated  into  a  set  of  rate  equations  which  have  the 
formZ : 

-Z-,j-  =  [N  N(Z-1 , j" )S(Z-1 , j" :Z,j)  +  N  N  (Z+I.j' ) 
dt  e  e 

[a (Z+1 ,j' :Z,  j)  ♦  Me8(Z+1,j' :Z,j)]  ♦ 

Z  N(Z, i)A(i , j)  +  N  Z  N(Z,i)X-1(i,j)  + 
i>j  i>j 

Ng  Z  N(  Z, i ) X(i ,  j )  ] 
i>j 

-N(Z,  j)[NeS(Z,j:Z-t.1,j'  )  +  Ng[a  (Z,  j  :Z-1 ,  j"  )  + 

N  8( Z , j : Z- 1 , j" ) ]  +  Z  A(j,i)  + 

6  i>j 

Ne  2  X_L(j,i)  +  Ne  Z  X(j,i)] 
i>j  i>j 

where  N(Z,j)  i3  the  population  density  of  ion  species  Z  in  level  j, 

S(Z-1,j";  Z,j)  is  the  collisional  ionization  rate  (cm  /sec)  from 

N(Z-1,j")  to  N(Z,J),  a  is  the  sum  of  the  radiative  and  dielectronic 

recombination  rates  (ar  +  a^r)  in  cm  /sec,  6  is  the  collisional 

6 

recombination  rate  (cm  /sec),  A(i,j)  is  the  spontaneous  decay  rate  from 

levels  i  to  J  (1/sec),  X(i,j)  is  the  collisional  excitation  rate  from  i 
3  -1 

to  j  (cm  /sec)  and  X  is  the  inverse  reaction,  collisional 
de-excitation . 

(Note  that  since  this  is  an  equilibrium  model,  the  left  hand  side  of  the 
above  equation  is  set  equal  to  zero.) 
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In  this  form  the  model  consists  of  an  infinite  set  of  coupled 
non-linear  equations,  one  for  each  available  state  of  the  ions  and  atoms 
present  in  the  plasma.  Fortunately,  there  are  some  methods  and 
procedures  available  for  truncating  the  number  of  states.  Though  we  are 
not  explicitly  considering  reduction  in  ionization  potentials  in  this 
model,  this  phenomenon  still  sets  an  upper  bound  on  the  number  of 


discrete  quantum  states  available,  i.e.,  we  only  need  model  those  states 
for  which  E(n)  <  E(Z)-AE(Z). 

In  our  model  we  have  considered  only  the  5-10  lowest  states  for 
each  ionization  stage  of  carbon,  including  neutral  carbon.  Although  a 
limited  number  of  states  have  been  considered  in  this  model,  they  should 
still  be  adequate  for  modelling  the  radiation  field,  specific  heats  and 
other  equations  of  state  information  needed  for  the  beam  -  target 
simulations. 

We  have  shown  that  the  CRE  model  is  applicable  over  the  range  of 

plasma  conditions  considered  in  this  work,  i.e.,  for  temperatures  between 

19  24  -3 

1  and  100  ev  and  ion  densities  between  10  and  10  cm  .  Both  the 
CRE  and  LTE  models  give  identical  results  for  collisionally  dominant 
plasmas.  Similarly,  the  CRE  model  exhibits  coronal  behavior  in  the  low 
density  limit.  This  ability  to  model  both  collisional  and  radiative 
plasmas  is  essential  for  a  comprehensive  ionization  dynamics  model. 

In  addition,  we  investigated  the  time  constraints  placed  upon  the 
overall  simulation.  In  particular  this  consists  of  a  consideration  of 
the  equilibrium  aspects  of  the  CRE  model  and  the  single  temperature  and 
fluid  assumptions  inherent  in  the  hydrodynamics  code.  In  all  cases,  we 
found  that  the  models  are  valid  for  the  timing  limitations  encountered 
in  this  study. 

We  have  calculated  temperature,  density  and  pressure  profiles  for 
evolving  targets.  We  find  that  the  deposition  region  heats  very  rapidly 
until  the  processes  of  electron  thermal  conduction,  radiative  loss  and 
beam  deposition  begin  to  equilibrate.  Once  this  occurs,  we  find  that 
the  pressure  across  the  width  of  the  deposition  region  remains  fairly 
constant  during  the  interaction.  The  pressure  is  responsible  for 
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compression  of  the  nondeposition  region  to  densities  as  high  as  five  to 
ten  times  solid  state.  We  provide  diagnostic  information  that  displays 
various  transition  energies  and  power  densities  during  the  ablation.  In 
addition,  the  partitioning  of  the  deposited  energy  into  internal, 
kinetic  and  radiative  components  is  examined.  This  model  can  be  applied 
to  a  variety  of  plasma  and  beam  conditions.  The  plasma  restrictions 
have  already  been  mentioned,  and  the  beam  conditions  must  be  such  that 
the  energy  and  power  density  is  sufficient  to  heat  the  target  to  at 
least  several  ev  (1  ev  represents  the  lower  bound  of  the  rate 
coefficient  table). 

An  examination  of  the  results  of  our  parameter  study  suggests  the 
direction  one  might  take  in  order  to  tailor  the  energy  partitioning  of 
the  radiation  field  or  ablation  region  for  a  specific  application.  For 
example,  if  carbon  Is  to  be  used  as  a  x-ray  laser  pump,  then  one  might 
want  to  optimize  the  radiation  conversion  efficiency  or  the  shape  of  the 
emission  profile.  In  our  work,  we  find  that  conversion  efficiencies  can 
be  as  large  as  90%.  We  also  find  that  the  thermal  emission  profile  is 
dependent  upon  both  the  energy  and  power  density  of  the  beam.  This  is 
illustrated  by  the  observation  that  the  hardness  of  the  radiation  field 
can  be  increased  by  either  using  less  energetic  protons  and  maintaining 
the  same  power  density  or  else  increasing  the  power  density.  The 
parameter  study  also  categorizes  the  dependence  of  internal,  kinetic  and 
radiative  energy  upon  beam  power  density  and  proton  energy.  The 
radiation  field  is  further  divided  into  relative  contributions  of  line, 
recombination  and  bremsstrahlung  radiation. 

We  also  investigated  the  effects  of  including  the  inner-shell  photon 
absorption  process.  This  is  the  only  radiation  transport  we  considered 
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PLASMA  "ENERGY" 
PARAMETERS 

INITIAL  ION 
ENERGY 

1  Mev 

INITIAL  ION 
ENERGY 

2  Mev 

INITIAL  ION 
ENERGY 

3  Mev 

ENERGIES  IN  UNITS  OF  { 103  J/CM2)  AND 

(  %  OF  TOTAL  BEAM  ENERGY)  ARE 

LISTED 

BEAM 

9.0  (100%) 

9.0  (100%) 

9.0  (100%) 

INTERNAL ( * ) 

1.9  (21%) 

2.7  (30%) 

3.2  (36%) 

KINETIC 

1.2  (13%) 

1.1  (12%) 

.95  (11%) 

RADIATIVE 

5.9  (65%) 

5.2  (58%) 

4.9  (54*) 

*** 

*** 

*** 

*** 

IONIZATION 

.87  (10%) 

1.0  (11%) 

1.3  (14%) 

LINE 

3.0  (32%) 

2.4  (27%) 

2.4  (27%) 

RECOMBINATION 

1.8  (20%) 

1.4  (15%) 

1.2  (14%) 

BREMSSTRAHLUNG 

1.2  (13%) 

1.4  (15%) 

1.3  (14%) 

LINE  SPECTRUM  ENERGIES  (%  OF  RADIATIVE  ENERGY) 

E<10ev 

3% 

7% 

13% 

10<E<20ev 

15% 

26% 

28% 

20<E<60ev 

31% 

13% 

8% 

60<E<100ev 

*** 

**  * 

*** 

100<E<200ev 

*** 

*** 

*** 

200<E<350ev 

*** 

*** 

*** 

E>350ev 

*** 

*** 

*** 

RECOMBINATION  SPECTRUM  ENERGIES  (%  OF  RADIATIVE  ENERGY) 

E<10ev 

*** 

1% 

5% 

10<E<20ev 

17% 

15% 

13% 

20<E<60ev 

13% 

9% 

6% 

60<E<100ev 

*•* 

*** 

*** 

100<E<200ev 

*** 

*** 

*  *  * 

200<E<350ev 

*** 

*** 

**  * 

E>350ev 

*** 

*** 

*** 

*  Ionization  energy  is  included  in  internal  energy 
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PLASMA  "ENERGY" 
PARAMETERS 

INITIAL  ION 
ENERGY 

1  Mev 

INITIAL  ION 
ENERGY 

2  Mev 

INITIAL  ION 
ENERGY 

3  Mev 

ENERGIES  IN  UNITS  OF  ( 105  J/CM2)  AND 

(  %  OF  TOTAL  BEAM  ENERGY)  ARE 

LISTED 

BEAM 

9.0  (100%) 

9.0  (100%) 

9.0  (100%) 

INTERNAL (* ) 

.89  (10%) 

.45  (  5%) 

.52  (  6%) 

KINETIC 

.49  (  5%) 

.47  (  5%) 

.45  (  5%) 

RADIATIVE 

7.7  (85%) 

8.2  (91%) 

8.1  (90%) 

*** 

*** 

*** 

*** 

IONIZATION 

.22  (  2%) 

.17  (  2%) 

.18  (  2%) 

LINE 

2.3  (25%) 

2.6  (29%) 

1.5  (17%) 

RECOMBINATION 

4.8  (53%) 

4.1  (45%) 

3.9  (43%) 

BREMSSTRAHLUNG 

.60  (  7%) 

1.5  (17%) 

2.6  (29%) 

LINE  SPECTRUM  ENERGIES  (%  OP  RADIATIVE  ENERGY) 


E<10ev 

*** 

*** 

*** 

10<E<20ev 

*** 

*** 

*** 

20<E<60ev 

2% 

6% 

11% 

60<E<100ev 

*** 

*** 

*** 

100<E<200ev 

*** 

*** 

*** 

200<E<350ev 

8% 

14% 

4% 

E>350ev 

14% 

11% 

3% 

RECOMBINATION  SPECTRUM  ENERGIES  (%  OF  RADIATIVE  ENERGY) 

E<10ev 

*** 

*** 

*** 

10<E<20ev 

*** 

*** 

*** 

20<E<60ev 

*** 

4% 

20% 

60<E<100ev 

8% 

20% 

19% 

100<E<200ev 

19% 

18% 

8% 

200<E<350ev 

5% 

*** 

**  * 

E>350ev 

30% 

7% 

1% 

Ionization  energy  is  included  in  internal  energy 


in  our  otherwise  "optically  thin"  model.  The  results  of  the  study  show 
that  as  much  as  10%  of  the  radiation  field  can  be  absorbed  by  the 
inner -shell  process.  Although,  for  an  optically  thick  application  this 
value  would  be  reduced  because  of  competition  with  valence 
photo-processes,  it  still  represents  an  important  absorption  mechanism 
for  recombination  radiation.  An  examination  of  the  temperature,  density 
and  pressure  profiles  show  that  inner-shell  photon  absorption  primarily 
affects  the  temperature  and  density  profiles  of  the  ablation  region. 

The  nondeposition  region  appears  to  be  insulated  from  its  effects; 
thermal  conduction  and  ion  shock  heating  remain  the  principal  heating 
mechanisms  in  this  region. 

Another  application  investigated  was  the  calculation  of  the  amount 
of  Ka  photon  energy  emitted  by  the  direct  interaction  of  the  beam 
protons  with  the  target  atoms.  We  found  that  under  optically  thin 
conditions,  the  energy  density  of  the  emission  is  several  orders  of 
magnitude  below  that  of  the  primary  radiation  field.  However,  since  the 
energy  of  the  Ka  is  larger  than  the  other  photons,  it  should  be  easily 
resolvable.  It  was  noted  that  this  situation  only  arises  when  the 
target  is  relatively  cold.  Otherwise,  there  will  be  a  significant 
component  of  high  energy  photons  created  by  valence  transitions  that 
will  also  be  absorbed  by  the  K  shell.  When  this  happens,  additional 
emission  is  produced  and  it  washes  out  the  beam  results.  Because  the 
spectrum  of  depends  upon  the  binding  energy  of  the  K  shell  electron 
of  the  emitting  atom,  it  represents  a  measure  of  the  ionization  state  of 
the  target. 

In  conclusion,  we  feel  that  this  model  contains  enough  essential 
physics  to  be  useful  for  analyzing  a  variety  of  processes  that  occur 


57 


during  the  ablation  of  a  carbon  target.  The  processes  which  we  examined 
in  this  investigation  included:  1)  the  effects  of  initial  beam 
conditions  on  plasma  parameters  such  as  kinetic,  internal  and  radiative 
energy,  2)  inner-shell  photon  absorption  and  3)  emission  due  to  the 
direct  interaction  of  the  beam  particles  with  target  atoms.  With  the 
implementation  of  the  suggestions  for  improvement  listed  above,  the 
model  should  be  adequate  for  providing  accurate  diagnostic  information 
for  a  variety  of  physical  problems. 
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II.  Experimental  Progress 

Ha)  Neutral / Ion  Interactions  with  a  Laser  Ablation  Plasma 

The  attenuation  of  beam  current  density  is  governed  by  the 
product  of  plasma  density  and  interaction  cross  section 
integrated  over  the  beam  particle  path  length  in  the  plasma, 
illustrated  in  Figure  1.  While  Coulomb  collisions  account  for 
the  attenuation  of  an  ion  beam,  neutral  f  articles  may  undergo 
charge  exchange,  momentum  scattering,  and  impact  ionization 
with  plasma  ions,  along  with  ionization  by  impact  with  plasma 
electrons.  The  Coulomb  cross  section  is  calculated  from  an 
expression  obtained  by  integration  of  the  Rutherford  scattering 
formula.  Of  the  neutral  interaction  cross  sections,  the  charge 
exchange  cross  section  appears  to  be  the  largest  and  is  to  date 
the  only  one  for  which  we  have  found  adequate  experimental  and 
theoretical  treatment  in  the  literature1  (H°  +  C+g 
interactions).  For  a  10  keV  hydrogen  atom  in  a  C+c*  plasma,  we 
estimate  for  the  total  neutral  interaction  cross  section 

^  ^  '  5  x  10~15  cm2 

cx 

Beam  attenuation  predictions  are  presented  in  Figure  2  for 
the  case  of  10  keV  protrons  or  hydrogen  atoms  traveling  20cm 
through  a  quiescent,  20eV  carbon  plasma  of  various  densities. 
Fokker-Planck  calculations  of  initial  beam  particle  slowing 
down  (1^)  and  diffusion  ( )  times  indicate  that  path  lengths 
much  greater  than  20  cm  would  be  required  for  the  beam 
particles  to  acquire  an  isotropic  velocity  distribution  and  to 
thermalize  with  the  plasma,  since  a  10  xeV  proton  will  travel 
20  cm  in  about  140  ns.  However,  since  the  Coulomb  cross 
section  is  so  large  for  the  conditions  considered,  all  beam 
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Ion  beam  attenuation: 


I(t)  =  I0  exp 


where  O’ =  Coulomb  scattering  cross  section 


Neutral  beam  attenuation: 


I(t)  =  I0  exp 


dl 


where  z  =  cross  section  for  charge  exchange/  <£x , 

+  momentum  scattering,  Os  , 

+  impact  ionization,  O^, 
with  plasma  ions  of  species  i, 
charge  z,  and  density  n.  ; 

<OeVe>  =  electron  impact  ionization  cross 

section  averaged  over  the  electron 
velocity  distribution ; 

Vb  =  beam  particle  velocity  ; 

lp  =  beam  path  length  in  the  plasma  . 


Figure  1  -  Plasma  attenuation  of  a  particle  beam. 


60 


ions  will  undergo  Coulomb  collisions  while  traversing  the 

plasma  and  hence  will  be  lost  from  the  beam.  With  10  keV 

hydrogen  atoms,  on  the  other  hand,  we  predict  significant 

attenuation  by  charge  exchange  interactions  to  occur  for  plasma 

12  -3 

line  densities  greater  than  about  10  cm 

With  the  assumption  of  cross  sections  which  are 

independent  of  space  and  time,  we  see  in  Figure  1  that  the  time 

behavior  of  the  beam  current  density  exiting  the  plasma  depends 

on  the  time  behavior  of  the  integral  of  plasma  density  over 

beam  particle  path  length  in  the  plasma,  which  is  defined 

in  terms  of  the  line  density  along  the  particle  beam  axis. 

An  expression  for  the  line  density  resulting  from  the 

production  of  a  laser  ablation  plasma  is  obtained  from  the 

2 

integration  of  Tallents'  expression  for  the  ablation  plasma 

ion  density,  presented  in  Figure  3.  The  laser  beam  is  taken 

to  be  perpendicular  to  both  the  particle  beam  axis  and  the 

target  surface.  The  buildup  and  decay  of  the  line  density 

after  firing  the  laser  is  also  presented  in  Figure  3. 

lib)  Duopigatron  Based  Experiments  Concerning  Neutral/Ion 
Beam  Interactions  with  Laser  Ablation  Plasmas 

During  the  past  year  our  experimental  efforts  have  been 

concentrated  toward  the  measurement  of  the  effects  of  an 

independently-produced  ablation  plasma  upon  the  penetration  of 

a  neutral  or  ion  beam  to  an  underlying  surface.  This  problem 

has  relevance  to  the  case  of  a  plasma  or  gas  cloud  surrounding 

an  outgassing  material  in  a  space  vacuum  environment. 


61 


•  10  keV  proton  beam,  C4  plasma* 

O’  ~  (Tc  (Coulomb  scattering),  KT  ~  20  eV 


•  1 0  keV  neutral  hydrogen  beam,  C4  plasma 


•Assumed  uniform  and  stationary  (  n  =  n  L  ,  L  =  20  cm  ) 

Figure  Beam  attenuation  predictions 


Ion  number  density  : 

(  G.  J.  Tallents,  Laser  and  Particle  Beams  (1983), 
V*  Vol.  1,  Part  1,  pp.  171-80. 


Line  density  along  particle  beam  axis,  located  a 
distance  za  from  target  surface: 


5  10  15  20  25  30  35  40 

Time  (jjs  ) 


Conditions:  z#  =  1  cm,  L  =  20  cm,  v°  =  5  X  106cm/s, 

V*  Z 

k  =  =  4,  Nt  =  1015  ions. 

r 

Figure 3  -  Self-similar  expansion  model  for  laser-produced 
target  plasma 
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Ilb.l)  Beaml ine  and  Vacuum  Systems 

The  duopigatron  based  experiment  has  been  completely 
disassembled  and  renovated  for  these  experiments.  What  was 
needed  for  these  interaction  experiments  was  a  highly 
collimated  ion  or  neutral  beam  to  interact  with  the  expanding 
carbon  ablation  plasma  produced  by  a  ruby  laser.  A  floor  plan 
showing  the  present  location  of  the  components  in  this 
experiment  is  shown  in  Figure  4.  Major  components  and  changes 
in  the  system  will  be  discussed  in  the  following  sections. 

The  neutral  beam  device  is  depicted  as  it  existed 
immediately  following  its  renovation  in  Figures  5  and  6.  Large 
port  gate  valves  are  used  to  isolate  components  of  the  system. 

The  beamline  expands  from  the  10.2  cm  diameter  opening  between 
the  neutralization  tank  and  the  vertical  gate  valve  to  the  15.2 
cm  diameter  opening  into  the  target  chamber,  which  is  a  29  cm 
diameter,  33  cm  long  stainless  steel  tank.  Following  vacuum 
leak  and  electrical  system  checks,  10  kV  hydrogen  beams  were 
extracted  and  monitored  to  verify  proper  operation. 

In  order  to  place  diagnostics  so  as  to  monitor  the  beam 
after  passage  through  a  target  plasma,  a  number  of  changes  were 
made  to  the  system.  First,  in  order  to  avoid  filling  the  target 
chamber  with  a  background  plasma  produced  by  passage  of  the  beam 
into  the  chamber,  and  to  better  collimate  the  particle  beam 
entering  the  chamber,  a  3  cm  aperture  was  placed  at  the  chamber 
entrance.  Better  results  were  obtained  by  adding  two  more 
apertures,  one  just  after  the  vertical  gate  valve  and  another 
just  after  the  second  pumping  station.  This  arrangement  provided 
bc':h  a  collimated  particle  beam  and  differential  pumping^ 
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Figure  4  -  Components  of  the  beam  interaction  experiment.  Shown  is  (1)  the  duopigatron. 

(2)  neutralization  cell,  (3)  bending  magnets,  (4)  target  chamber,  (5)  particle  detector,  (6)  isolation 
transformer,  (7)  isolated  power  supplies,  (8)  controls,  (9)  ruby  laser,  (10)  HeNe  laser,  (11)  gas  bottles 
Not  shown  is  the  hioh  voltaae  power  supplu  located  in  a  separate  interlocked  room. 
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Figure  5  -  Schematic  illustration  of  duoPIGatron  type  positive 
ion  source  with  accelerating  grids  and  neutralizing 
gas  cell.  66 
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Duopigatron  and  beamline  as  constructed  at  the  beginning  of  this  reporting  period. 


Pressure  in  the  target  chamber  remains  below  10-^Torr  while  the 

_3 

pressure  in  the  duopigatron  end  rises  to  nearly  10  Torr  when  a 
beam  is  extracted.  Second,  since  the  experimental  targets  are  to 
be  located  in  the  plane  defined  by  the  11  ports  arranged 
circumferentially  near  the  end  of  the  chamber,  a  10  cm  I.  D. 

Pyrex  cross  was  attached  to  the  target  chamber  end  flange  to 
allow  the  mounting  of  diagnostics  downstream  of  the  target. 
Initial  experiments  using  thi  arrangement  indicated  the  need 
to  move  the  diagnostics  further  away  from  the  target,  due  to 
target  plasma  streaming  and  particle  beam  timing,  and  so  a 
second  Pyrex  section  was  placed  between  the  first  one  and  the 
chamber  end  flange.  Further  studies  indicated  the  need  for 
additional  apertures  to  be  placed  between  the  target  plasmc  and 
the  Pyrex  system  along  with  a  reduction  of  the  size  of  the 
aperture  at  the  entrance  to  the  target  chamber .  The  aperture 
into  the  Pyrex  system  made  it  necessary  to  install  a  third 
vacuum  pumping  station.  The  target  end  of  the  neutral  beam 
device,  as  of  the  end  of  this  reporting  period,  is  illustrated 
in  Figure  7. 

IIb.2)  Ruby  Laser  System 

A  target  plasma  is  produced  in  the  target  chamber  by 
focusing  the  output  of  a  Korad  K-15  ruby  laser  onto  a  solid 
target  mounted  within  the  chamber  on  the  end  of  a  rotatable 
feedthrough.  A  29  cm  focal  length  lens  is  mounted  external  to 
the  vacuum  system  and  the  laser  beam  passes  through  a  glass 
window.  A  prism  is  used  to  turn  the  laser  beam  and  a  HeNe 
laser  is  used  to  align  the  laser  system  components,  which  are 


SI— H 


arget  end  of  beamline,  (dimensions  in  cm) 


shown  in  Figure  4.  A  pockels  cell  allows  conversion  of  the 
laser  output  from  the  relaxation  mode  to  the  Q-switched  mode. 


IIb.3)  Diagnostics 

Beam  attenuation  data  is  obtained  by  the  use  of  a  Faraday 
cup  type  of  particle  detector,  illustrated  in  Figure  8.  In 
this  device,  an  energetic  particle  (ion  or  neutral)  strikes  the 
collection  plate,  which  then  emits  secondary  electrons.  The 
secondary  electrons  are  collected  by  the  innermost  cup,  which 
is  typically  biased  at  +15  VDC.  The  outer  cups  define  the 
collector's  acceptance  angle  and  help  to  shield  the  innerccup 
from  any  charged  background  particles.  Ions  collected  on  the 
emitting  plate  and  electrons  leaving  the  plate  produce  a 
measurable  voltage  signal  across  the  resistor  shown  in  the 
figure . 

Other  beam  parameters  monitored  during  each  shot  are:  1) 
Arc  current,  I  ,  the  current  flowing  from  cathode  to  anode  in 
the  duopigatron;  2)  Drain  current,  Idrajn»  the  current  drawn  by 
the  accelerator  grid;  and  3)  Accelerator  voltage,  V  .,  the 
voltage  applied  to  the  accelerator  grid.  Arc  current  is 
obtained  by  measuring  the  voltage  across  a  shunt  resistor, 
drain  current  is  obtained  through  the  use  of  a  Pearson  coil  on 
the  transmission  line  between  the  high  voltage  capacitor  and 
the  accelerator  grid,  and  the  voltage  on  the  accelerator  grid 
is  measured  with  a  high  voltage  probe. 


IIb.4)  Ion  and  Neutral  Beam  Probe  Measurements 


particle  cetector  as  mounted 


experiment.  In  preparation  for  the  target  interaction 
experiments,  a  Pyrex  tee  was  attached  to  the  target  chamber  end 
flange  and  the  detector  was  mounted  on  a  movable  feedthrough. 

A  ruled  marker  was  used  to  indicate  the  cup  position  relative 
to  the  teamline  axis  and  the  cup  was  repositioned  between  shots 
to  scan  across  the  beam  diameter. 

Typical  oscilloscope  traces  showing  the  timing  and  beam 
parameters  for  20  amp  arc  discharges  and  an  accelerating 
voltage  of  10  kV  are  presented  in  Figure  9. 

Figure  10  presents  the  results  of  a  vertical  scan  of  beams 
after  having  passed  through  the  single  3  cm  aperture  for  both  0 
and  5  amps  of  bending  magnet  current  (0  and  approx.  180  Gauss). 
The  factor  of  three  decrease  in  the  peak  signals  for  such  a  low 
bending  magnet  current  indicates  that  the  beam  ion  fraction  is 
quite  high  and  the  degree  of  beam  divergence  is  indicated  by 
the  width  of  the  curves:  For  these  non-optimal  conditions, 
more  than  70%  of  the  beam  particles  entering  the  target 
chamber  are  ions  and  the  beam  expands  from  3  cm  to  about  6 
cm  in  diameter  in  a  distance  of  37  cm.  Fixing  the  cup 
position  at  2  cm  from  the  beamline  axis  and  varying  the 
current  through  the  bending  magnets  yielded  the  results 
presented  in  Figure  11.  This  shows  that  20  amps  (approx. 

700  Gauss  on  axis)  removes  all  but  the  heaviest  of  beam 
contaminants  and,  again,  that  only  about  30%  of  the  beam 
particles  are  charge  neutral. 

Figures  10  and  11  were  obtained  from  operation  of  the 
source  with  a  relatively  high  hydrogen  gas  feed  rate  (gas 
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Figure  9  -  Typical  hydrogen  beam  parameters 
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valve  setting  =  100).  In  an  effort  to  increase  the  beam 
neutral  fraction  and  to  improve  beam  collimation,  two 
additional  3  cm  diameter  apertures  were  placed  in  the 
beamline  as  shown  in  Figure  7.  With  the  detector  placed  on 
the  beamline  axis,  10  kV  hydrogen  beams  were  again  extracted 
from  20  amp  arc  discharges,  this  time  varying  the  gas  feed 
rate.  Results  for  bending  magnet  currents  of  0  and  20  amps 
are  presented  in  Figure  12,  which  show  that  the  neutral 
fraction  of  the  beam  increases  as  the  gas  feed  rate 
decreases,  as  expected,  since  there  are  fewer  background 
particles  available  for  the  neutrals  to  charge  exchange  with 
(back  into  ions)  in  the  beamline.  However,  there  is  a  lower 
bound  on  the  gas  feed  rate  and  gas  pulse  length,  below  which 
there  is  no  source  plasma  buildup.  This  limit  appears  to  be 
at  a  gas  valve  setting  of  approximately  50  and  a  gas  pulse 
length  of  about  0.3  seconds  for  the  aperture  placement  shown 
in  Figure  7.  Discharges  at  this  operating  point  are  very 
unstable. 

Finally,  with  the  gas  valve  setting  at  60  and  gas  pulse 
length  of  0.3  seconds,  the  cup  was  again  scanned  across  the 
beam  diameter  to  obtain  the  results  shown  in  Figure  13. 
Although  the  apertures  have  reduced  the  beam  current  entering 
the  target  chamber,  these  results  clearly  show  that  neutral 
fraction  under  these  conditions  is  now  on  the  order  of  90%  and 
that  the  beam  collimation  is  somewhat  improved. 
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IIb.5)  Ablation  Plasma  Production  and  Diagnostics 

The  production  of  an  ablation  plasma  involved  focussing 
the  ruby  laser  output  onto  the  surface  of  a  solid  target:  A 
plume  of  ionized  target  material  expands  from  the  surface  at  a 
velocity  on  the  order  of  10^  cm/sec.  The  total  number  of 
particles  emitted  and  their  energy  depends  on  the  energy 
delivered  by  the  laser.  Initial  use  of  the  ruby  laser  involved 
its  operation  in  the  relaxation  mode,  which  can  deliver  up  to 
about  1.2  Joules.  However,  in  order  to  facilitate  the 
correlation  of  laser  output  with  diagnostic  signals,  the  laser 
has  been  operated  during  this  period  in  the  Q-switched  mode, 
which  delivers  up  to  about  0.3  Joule  with  a  pulse  width  of  40- 
60  ns.  Oscillograph  traces  from  PIN  diode  detection  of  the 
laser  output  in  both  modes  are  presented  in  Figure  14. 

Prior  to  this  reporting  period,  carbon  plasmas  were 
produced  using  the  relaxation  mode  output  and  attempts  were 
made  to  characterize  them.  These  initial  studies  involved  the 
use  of  flat  copper  plates,  biased  at  various  voltage,  to 
collect  the  charged  plasma  particles.  With  the  plates  located 
1.55  and  2.55  cm  from  the  face  of  the  graphite  target,  it  was 
determined  that  the  electron  temperature  is  in  the  range  16  < 
kT(eV)  <  24,  that  the  maximum  ion  velocity  is  about  2  x 
lO^cm/s,  which,  for  carbon  ions,  yields  an  average  expansion 
energy  of  6.2  eV. 

More  recently,  Q-switched  laser  produced  carbon  plasma  has 
been  studied  through  the  use  of  a  Langmuir  probe  in  order  to 
obtain  an  estimate  of  the  plasma  electron  temperature.  With 
the  laser  beam  axis  perpendicular  to  the  graphite  target,  the 


probe  was  placed  20  degrees  off  axis  and  4.3  cm  from  the  focal 
spot.  The  results  of  these  experiments  are  presented  in  Figure 


15  and  calculations  yield  electron  temperatures  of  8.7,  11.1, 
and  10.7  eV  at  2,  5,  and  10  us,  respectively. 

We  have  also  looked  at  the  ablation  plasma  emission 
spectra  by  replacing  the  Pyrex  vacuum  chamber  with  a  large 
glass  window  on  the  target  chambe  and  flange  and  focusing  the 
light  onto  the  entrance  slit  of  an  optical  multichannel 
analyzer  (OMCA) .  Spectra  obtained  from  three  separate  shots 
are  presented  in  Figure  16.  Although  we  were  unable  to 
calibrate  for  intensities  to  obtain  line  ratios  and  obtain  plasma 
temperature,  these  spectra  show  that  we  are  producing  carbon  ions 
within  a  charge  state  of  at  least  +2. 

IIb.6)  Beam/Plasma  Interaction  Experiments 

We  have  measured  changes  occurring  when  energetic 
hydrogen  particles  interact  with  a  carbon  ablation  plasma. 

Timing  is  such  that  the  laser  is  Q-switched  after  the 
particle  beam  is  well  established.  Our  initial  experiments 
exhibited  a  low  signal  to  noise  ratio.  Further 
investigations  identified  more  than  one  source  for  the 
interference  received  by  the  detector  signal.  First,  it  was 
found  that  the  high  frequency  signal  occurring  when  the 
laser  is  Q-switched  was  due  to  the  detector  circuit  acting 
as  an  antenna,  receiving  the  noise  generated  by  the  pockels 
cell  circuit.  Second,  it  was  found  that  the  EMI 
(electromagnetic  interference)  problem  was  aggravated  by 
ground  loops  in  the  system.  Third,  it  was  found  that  the 
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detector  was  too  close  to  the  plasma,  that  there  was 
significant  streaming  of  plasma  into  the  detector. 

Adding  a  second  Pyrex  piece  to  place  the  detector 
further  from  the  target  chamber  helped  to  reduce  the  third 
source  but  not  to  eliminate  it,  and  so  small  apertures  were 
placed  between  the  target  chamber  and  the  Pyrex  system, 
requiring  the  addition  of  the  third  pumping  system.  Next, 
ground  loops  were  eliminated  and  attempts  were  made  to 
shield  from  the  EMI  by  placing  copper  sheet  and  RF  shielding 
tape  around  the  pockels  cell,  which  significantly  reduced 
the  magnitude  and  duration  of  the  noise  picked  up  by  the 
detector  circuit. 

With  the  noise  levels  reduced  to  the  point  where  we 
could  make  meaningful  measurements  of  the  small  change  in 
detector  current  following  the  firing  of  the  laser,  the  full 
beam/plasma  interaction  experiment  was  conducted  using  10 
keV  hydrogen  beams  and  Q-switched  laser  pulses.  Typical 
attenuation  data  for  the  case  of  no  bending  magnet  current 
(neutral  +  ion  beam)  and  for  18  amps  of  bending  magnet 
current  (neutrals  only)  are  presented  in  Figure  17.  It  is 
seen  that  as  much  as  12%  of  the  beam  current  is  lost  due  to 
beam  interactions  with  the  plasma. 

Data  from  a  "neutrals  only"  shot  is  presented  in  Figure 
18  along  with  a  best  fit  of  the  self-similar  expansion 
model.  The  model  yields  a  value  of  the  maximum  line  density 
which,  for  a  maximum  attenuation  of  approximately  13%, 
agrees  quite  well  with  the  simple  predictions  given  in 
Figure  2. 
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Figure  1 7  -  Typical  attenuation  data. 
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Model  :  Self-similar  expansion  normalized  to  the  data  at 
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Best  fit  obtained  with  k  =  4,  v*  =  3  x  106  cm/s 


Figure  1 8  -  Neutral  beam  attenuation,  experiment  vs  mode). 


Current  work  includes  providing  the  means  by  which  we 
can  discriminate  between  beam  ions  and  neutrals  in  order  to 
correlate  shots  such  as  are  shown  in  Figure  17.  We  are  also 
investigating  the  discrepancy  between  data  and  model  in  the 
first  5  us  after  firing  the  laser,  seen  in  Figure  18.  The 
ion  fraction  with  the  magnetic  field  on  is  less  than  1%. 

For  the  case  of  no  bending  magnet  the  ion/neutral  fraction 
is  estimated  to  be  13%:87%. 

He)  Fabrication  of  an  Applied  -  B  Ion  Diode  for  the  Febetron 
In  order  to  achieve  higher  energy  and  higher  current 
ion  beams  for  our  interaction  experiments  we  have  fabricated 
an  applied  B  (radial)  magnetically  insulated  ion  diode  which 
is  compatable  with  our  existing  Febetron  Pulser.  The  higher 
energy  beam  will  permit  us  to  obtain  higher  energy  scaling 
data  compared  to  the  20  kV  duopigatron  data,  whereas  the 
higher  beam  current  should  improve  the  signal  to  noise  ratio 
of  our  particle  detectors.  The  Febetron  parameters  are: 

Voltage  Maximum  =  500  kV 

Current  Maximum  =  1.4  kV 

Pulselength  =  400  nsec,  full  width  on  current 

A  schematic  illustration  of  the  applied  -  B  diode  as 
configured  for  the  Febetron  is  depicted  in  Figure  19.  (This 
design  is  similar  to  those  described  in  references  4  through 
9.)  The  Febetron  will  initially  be  run  in  positive  polarity 
with  a  flashboard  anode  and  a  carbon  brush  cathode.  Figures 
20  and  21  give  more  detailed  views  of  the  recently 


constructed  ion  diode  magnet  coils  and  vane  structure.  The 
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flashboard  anode  (Fig.  22)  employs  a  lucite  plate  with 
imbedded  copper  pins. 

An  important  consideration  in  the  design  of  this  diode  was 
the  necessity  of  achieving  a  radial  magnetic  field  which 
exceeds  the  critical  insulating  field.  The  critical  insulating 
field  for  a  1-D  planar  diode  can  be  calculated  from  the 
expressions  given  in  references  10-13.  Assuming  that  we  wish 
to  operate  at  1.5  times  the  critical  field  the  following  table 
gives  the  required  performance  of  our  ion  diode: 

V(kV)  B  (Gauss)  1.5  B  (Gauss) 

crit  crit 


50 

772 

1158 

100 

111? 

1677 

200 

1650 

2475 

300 

2101 

3152 

400 

2516 

3775 

500 

2911 

4366 

The  magnitude  of  these  fields  dictated  that  the  magnet 
coils  needed  to  be  operated  in  the  pulsed  mode.  The 
capacitive  discharge  magnet  power  supply  depicted  in  Figure 
23  was  constructed  for  this  purpose.  Details  of  the 
ignitron  firing  circuit  are  given  in  Figure  24.  Typical 
data  for  the  magnet  current  trace  and  integrated  B  dot  loop 
signals  are  given  in  Figure  25.  The  measured  radial 
magnetic  field  (Figs.  26  and  27)  easily  exceeds  the  critical 
insulating  field  for  the  magnet  and  power  supply  which  we 
have  constructed. 
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Figure  23  -Pulsed  Magnet  Supply 
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Ftg.25-  Top  trace:  example  of  self-integrated  B  dot  probe  (r  =  center 
of  gap,  z  =  3.5  cm)  100  turns,  A  =  4.95  x  10-5m2,  RC  =  22  ms,  V  =  2  KV 
Note  that  the  trace 'for  B(t)  follows  the  current  very  well  signifying 
good  integration  of  the  signal. 

Bottom  trace:  Magnet  current  measured  through  Pearson  coil  0.05  v/amp 
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Figure  27 
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